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Abstract 

The first radiative correction to the Casimir energy of a perfectly con- 
ducting spherical shell is calculated. The calculation is performed in the 
framework of covariant perturbation theory with the boundary conditions 
implemented as constraints. The formalism is briefly reviewed and its use 
is explained by deriving the known results for two parallel planes. 

The ultraviolet divergencies are shown to have the same structure as 
those for a massive field in zeroth order of a. In the zeta-functional regu- 
larization employed by us no divergencies appear. 

If the radius of the sphere is large compared to the Compton wavelength 
of the electron the radiative correction is of order a/(R 2 m e ) and contains 
a logarithmic dependence on m e R. It has the opposite sign but the same 
order of magnitude as in the case of two parallel planes. 

1 Introduction 



The Casimir effect is one of the basic effects in Quantum Electrodynamics (QED) 
and provided the first example of Nullpunktsenergie in a field theory. Proposed in 
1948 and qualitatively verified in 1958, it has only recently been tested quantita- 
tively |]]. Based on the attractive force that is exerted on the conducting plates, 
Casimir even proposed a model for the electron. This model looked quite promis- 
ing until Boyer showed a repulsive force in the first calculation for a sphere in 
1968. In a large number of calculations for different geometries the Casimir force 
was found to be repulsive for configurations whose sizes in different directions are 
close to one another (a cube, for instance) and attractive when at least one of the 
sizes is much larger than the others (a long cylinder, for instance). Nevertheless, 
a general understanding of this property is still missing. 
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Radiative corrections have rarely been considered because the expected effect, 
being proportional to the fine structure constant a and at least to the ratio of 
the Compton wavelength A c of the electron to a typical geometric size L of the 
boundaries, is too small to be directly observable. Nevertheless, there is a general 
interest in their consideration, mainly having in mind the applicability of the 
general methods of quantum field theory For example, its covariant formulation 
in the presence of explicitly non-covariant boundaries and the generalization of 
the perturbation expansion of a gauge theory are among the important issues 
to be addressed in this context. Boundary conditions necessitate a modification 
of the renormalization procedure as well. The interaction of the quantum fields 
with the boundary leads to ultraviolet divergencies in vacuum graphs that cannot 
simply be discarded by normal ordering. Last but not least the question for the 
order of the geometrical effects is raised, i.e. the question for the leading power 
of A c / L contributing to the effective action. Additional attention to radiative 
corrections arises in the bag model of QCD where they are not a priori negligible. 

The calculation of the ground state energy can be performed starting from 
the basic relation 

(n) 

where Ei n \ are the one particle energies of the considered system and the sum 
runs over the corresponding spectrum. This form, modified by some regulariza- 
tion (e.g., the zeta functional one with E r Q cg = (1/2) J2( n ) E \n) S i ( s > 3 / 2 )), is 
best suited for most of the calculations for different boundary conditions and ge- 
ometries. /,From the point of view of general quantum field theory it is connected 
with the vacuum expectation value of the energy- momentum tensor < 0|T MZ ,|0 > 
by Eq = J dx < 1 Too 1 > (leaving aside modifications due to translational invari- 
ant directions). By means of the well known relations T 00 = -£-(p(x)-£-(f(x) (real 

scalar field), < O\0{x)0(y)\O >= iD(x,y) and < 0|f 00 |0 >= iJLJL-D(x,y)\ v=x it 
is qualified as a one loop contribution and can be represented as the Feynmann 

graph \^^y' w ^ ere ^ * s ^ ne ver ^ ex factor corresponding to T^ v . Using this 
language it is straightforward to incorporate radiative corrections as higher loops. 
For instance, in <^ 4 -theory we would have <^^^(^)- I n QED there are two contri- 
butions to the energy-momentum tensor, one resulting from the electromagnetic 

field w ) and one from the spinor field <S> )• The corresponding two 
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loop contributions are (V) k t and 

Currently there is no established formalism for the calculation of radiative 
corrections in the presence of boundaries, although a number of attempts have 
been undertaken || [|. In the present paper we use the formalism of general 
covariant perturbation theory, including the standard theory of renormalization 
with the 'minimal' modifications needed and calculate the first radiative correc- 
tion to the Casimir energy for a perfectly conducting spherical shell. The basic 
ideas are taken from an earlier work ||, however the time passed shows that the 
required modifications are not simple. 

The conductor boundary conditions the electromagnetic field has to fulfil on 
some boundary S read 

n"F; u (x) Us = , (2) 

where F* = e^ vol pF af5 is the dual field strength tensor and n M is the (outer) 
normal of S. Being the idealization of a physical interaction (with the conductor), 
they are formulated in terms of the field strengths and thus gauge-invariant. 
Now, for well known reasons, it is desirable to perform the quantization of QED 
in terms of the gauge potentials A^{x). Obviously, the boundary conditions (§) 
do not unambiguously imply boundary conditions for all components of A^(x) as 
it is required in order to obtain a selfadjoint wave operator. 

There are two possibilities to proceed. The first one is to impose boundary 
conditions on the potentials in such a way that the conditions (Q) are satisfied 
and a selfadjoint wave operator is provided. Electric and magnetic boundary 
conditions are commonly used for this purpose. These conditions are stronger 
than @ but they are not gauge invariant. When now requiring BRST invariance 
the ghosts become boundary dependent too and contribute to physical quantities 
like the Nullpunktsenergie. This was first observed in the papers || and has 
later also been noticed in connection with some models in quantum cosmology 
[0. The common understanding is that the ghost contributions cancel those from 
the unphysical photon polarizations, for recent discussions see || [| . 

In a second approach one considers the boundary conditions (|2]) as constraints 
when quantizing the potentials A^(x) as it was first done in |J. In that case, 
explicit gauge invariance is kept. There is no need to impose any additional con- 
ditions. The conditions (§) appear to be incorporated in a 'minimal' manner. In 
that respect, this second approach resembles the so called dyadic formalism which 
has succesfully been used in the calculation of the Casimir energy in spherical 
geometry [|1(] . 

There are two ways to put this approach into practice. The first way is to 
solve the constraints explicitly. For this purpose one has to introduce a basis of 
polarization vectors E s (instead of the commonly used e s A such that only two 
amplitudes (in our notation those with s — 1, 2) of the corresponding decompo- 
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sition 



a. q (x 



(3) 



of the electromagnetic potential have to satisfy boundary conditions. The other 
two amplitudes (in our notation those with s = 0, 3) remain free. In the case 
where the surface S consists of two parallel planes such a polarization basis can 
be constructed explicitly |TT ] 
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(E^g^El = g st ). Inserting the decomposition (|3|) with these polarization vectors 
into the boundary conditions (@), we find 



ai(x) = 0,2(0:) = for x G S , 



(5) 



whereas the amplitudes a s (x) with s = 0, 3 are unaffected by the boundary 
conditions. In this way, roughly speaking half of the photon polarizations feel the 
boundary (in ]T2] they have been shown to be the physical ones in the sense of the 
Gupta-Bleuler quantization procedure) and half of them do not. However, such 
an explicit decomposition, which simultaneously diagonalizes both, the action and 
the boundary conditions can only be found in the simplest case. The problem is 
that for a non-planar surface S the polarizations E^ become position dependent. 
To the authors' knowledge, their existence in the general case is still an open 
issue. 

Here we follow the second way to realize the boundary conditions as con- 
straints, namely we use Lagrange multipliers or, equivalently, restrict the inte- 
gration space in the functional integral approach. The generating functional of 
the Green functions in QED reads 



Z(J, rj, V ) = C j VA Vip Vip Y\.n$ (^F; u (x)) exp iS V, i> 
J v xes 

+ i dx \An(x) J^(x) + i/j(x)rj(x) + rj(x)tp(x) 



(6) 



where the integration runs over all fields with the usual asymptotic behavior. 
The functional delta function restricts the integration space to such potentials 
Afj, that the corresponding field strengths satisfy the boundary conditions (0). 
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This approach had been used in [EJ in the case of plane parallel boundaries. It 
was shown to result in a new photon propagator and an otherwise unaltered 
covariant perturbation theory of QED. As it will be seen in the next section, the 
boundary conditions (Q) appear to be incorporated with a 'minimal' disturbance 
of the standard formalism, completely preserving gauge invariance (as well as the 
gauge fixing procedure) and Lorentz covariance as far as possible. An additional 
advantage of this approach is the fact that the resulting formulas for practical 
calculations turn out to be rather simple. 

The spinor field deserves a special discussion with respect to the boundary 
conditions. We do not impose boundary conditions on the electron and consider 
the electromagnetic field and the spinor field on the entire Minkowski space with 
the conducting surface placed in it. In the case of S being a sphere we thus 
consider the interior and the exterior region together. In general, the surface S 
need not be closed. Only the electromagnetic field obeys boundary conditions on 
the surface S. The electron penetrates it freely, it does not feel the surface. As 
for a physical model one can think of a very thin metallic surface which does not 
scatter the electrons but reflects the electromagnetic waves. If the thickness of 
the metallic surface (e.g. 1 /im) is small compared to the radiation length of the 
metal (e.g. 1.43 cm for copper |18[), this approximation is well justified. However, 
the radiative corrections are of order (a\ c /L) with respect to the Casimir force 
itself and thus too small to be directly observable. A discussion of the validity of 
this approximation seems therefore a bit academic. 

We remark, that the situation is different in the case of the bag model. The 
boundary conditions of the gluon field and of the spinor field are connected by 
means of the equation of motion (this is because the field strength tensor enters 
the boundary conditions rather than the dual field strength tensor in (^])). Also, 
we cannot a priori expect the radiative corrections to be very small. However, 
a detailed discussion of this subject goes beyond the scope of the present paper. 
Nevertheless, we would like to regard the investigation of the radiative corrections 
in QED as a necessary step towards the corresponding calculations in QCD. 

For the needs of renormalization we consider the quantum fields embedded in 
a classical system, given by the surface S with a corresponding classical energy. 
For example, the classical energy of a sphere takes the form 

£ class = pV + aS + FR + k + - , (7) 

ti 

where V = fvri? 3 and S = 4nR 2 are the volume and the surface area, respectively. 
The basic idea is to renormalize the divergencies of the ground state energy by a 
redefinition of the corresponding parameters p, a, F, k, h in the classical part of 



the system. For recent comprehensive discussions we refer to the papers [13, 14]. 
The general structure of the potentially divergent contributions for an arbitrary 
surface S is completely known from the heat kernel expansion. Explicit formulas 



for the sphere are given in, e.g., JT5]. When one considers the quantum fields 



5 



in the interior as well as in the exterior of the sphere (as we do), the divergent 
contributions with odd powers of the radius R cancel. Therefore, it is sufficient 
to keep only the two contributions 

^class = aS + k (8) 

in the classical energy. Hence, as a hypothesis we expect divergent contributions 
proportional to R 2 and R° (i.e., proportional to a constant) when including the 
radiative correction. Indeed, this turns out to be the case (see section 5). The 
corresponding redefinition of a and k now depends on the fine structure con- 
stant a and the electron mass m (see eqs. ([f4]) and (|75|) below). Besides these 
boundary induced divergencies, the ordinary free Minkowski space vacuum con- 
tributions which are usually discarded by normal ordering, are expected to appear 
as well. In line with the common interpretation, their removal is understood as 
a renormalization of the cosmological constant. 

It is interesting to observe the following fact. For dimensional reasons, the 
divergent contributions @ to the ground state energy are proportional to positive 
powers of the mass, a ~ m 3 and k ~ m. Consequently, they vanish for massless 
fields, e.g. the pure electromagnetic field. However, this well known picture is not 
preserved when including the first radiative correction. Due to the interaction a 
dimensional parameter, namely the electron mass, is introduced and the divergent 
contributions proportional to R 2 and R° do no longer vanish. The necessity to 
consider the general strucure (^) of the classical energy in QED becomes evident. 

In the case of plane parallel plates of separation L, the ground state energy 
Eq including the radiative correction Eq to leading order in A c / L is given by [|| 

^plates = ^(0) F {1) _ VT 2 1 7T 2 «A C 

E ° ~ E ° +E ° ~ 720L3- + 2560 ^ 

Let us remark that the first loop correction is of order a and that the first 
possible contribution in the small ratio X c /L is present. In fact, once this is 
known, the remaining task consists in the calculation of the number 7r 2 /2560. 
The calculations that have appeared in the literature, show different results. Xua 
implements periodic boundary conditions for the potentials. In his result, all 
powers of A c /L are absent leaving only an exponentially small contribution. In 
the paper by Kong and Ravndal |2| the same boundary conditions as here are 
studied with the help of a different method, namely some effective field theory 
and the Euler-Heisenberg Lagrangian. Their result, E^ ~ (a 2 /L 3 )(A c /L) 4 , is of 
one order a smaller than (H). The renormalization has also been discussed in the 
paper |]. 

In view of the result (§) the radiative correction in leading order of A c / R for 
a conducting sphere of radius R could be expected as 

^sphere 0.092353 a\ c 
= — 2R — C ' ~R2" ' ( 10 ) 
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where the first term is Boyer's result (in the formulation of fllQfl). The correction 
is calculated in section 5. We find 

= -7.5788 ■ 1(T 4 ^ \n(mR) - 6.4833 ■ 1(T 3 (11) 
R l R 2 

Besides the term proportional to 1/R 2 it exhibits an additional logarithmic de- 
pendence ln(mR)/R 2 (m is the mass of the electron). 

The result (Q) allows for an interesting intuitive interpretation. Since no 
boundary conditions have been imposed on the spinors, a photon can pass the 
surface via the creation and subsequent annihilation of a virtual electron-positron 
pair. So the separation of the plates is effectively enlarged due to the radiative 
corrections. In this sense, the result ([5]) can be viewed as a 'renormalization' of 
the distance L separating the plates 

to leading order in X c /L. This was first noticed in As it was shown in 



11], the same 'distance renormalization' results when considering the radiative 
corrections to the photon states between the plates. This picture can also be 
adopted in the case of a conducting sphere. The radiative correction increases 
the effective radius R according to 



R^R 



R ^ 



(13) 



(to leading order in X c /R). 

The paper is organized as follows. In the next section we present the pro- 
cedure of quantization with constraints and in the third section we derive the 
general structure of the radiative corrections. In the fourth section we rederive 
the results for parallel planes in order to explain the use of the formalism in 
details. The fifth section contains the main part of the paper, namely the calcu- 
lation of the radiative correction for spherical geometry. We conclude the paper 
with a summary and two technical appendices. 



2 Quantization 



We start the quantization procedure with the generating functional Z(J,rj,r)) 
represented by the path integral (H). The perturbation expansion is obtained by 
means of 

'J S_ _Sj 

SiJ 1 Sir]' 5ifj 



Z(J,r],r}) = exp 



iS\r 



(14) 



with Si nt (A, ip, ijj) = e J dx ip(x)A(x)ip(x) and the problem is reduced to the 
calculation of the generating functional for the non-interacting theory. 
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Before proceeding further we rewrite the boundary conditions (|2|) in the fol- 



lowing way. Let E s [ 
properties 



x 



(s = 1, 2) be the two polarization vectors in (]3|) with the 



d 

dx„ ^ v 
n»E;{x) 








(15) 
(16) 



for x G S (assuming d^n u (x) = d u n^{x)). They span a space of transversal 
vectors tangential to the surface S. Note that due to (|16|), there is no derivative 
acting outside the tangential space, i.e. no normal derivative in ([15]). Without 
loss of generality, we assume the normalization EfJ g^E^ = —5 st . We remark 



that the transformations — > + d tl ip(x) 



ft 

and A, 



A^ + n^ip(x) respect 



the boundary conditions (@), i.e. if A^(x) satisfies the boundary conditions, the 
transformed potential also does. The invariance under the first transformation 
simply says that the boundary conditions are gauge independent. The second 
means that the projection of A^ onto the normal of S is unaffected by the 
boundary conditions. We therefore conclude that the boundary conditions (0) 
can equivalently be expressed as 



EiA»(x) Us = (a = 1,2). 



(17) 



Explicit examples for such polarization vectors E s (besides E s with s — 1,2 in 
(HI) for parallel plates) are 
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El 



,Q2 + Q2 



( id X3 \ 



V id xo J 

t d X0 L \ 

R\ U X U X 3 J 
R\ U X U X 3 ) 

V 9 X3 L J 



-dl+dlJ-d 2 X0 +d 2 xs +L 2 



for a cylinder with axis along the a^-direction, radius R and the orbital momen- 
tum operator L = —x 2 d Xl + X\d X21 and 



El 



E, 




iRd X0 (n xL) I JJ?J-R?dl-L 2 



(19) 



for a sphere of radius R with normal vector n = x/\x\ and angular momentum 
operator L = ix x d x . The properties (|TB|) and (|TjJ) as well as the normalization 



S 



are easily checked. Evidently, these explicit polarization vectors do not contain 
normal derivatives. It should be noticed that the polarization E x „ in the examples 
(|18|) and (|19|) coincides with that of the standard approach (and the dyadic 
approach) whereas E 2 contains the combination n x L instead of V x L. 

Now, with the boundary conditions in the form ( |T7D we rewrite the delta 
functions in (§) and represent them as functional Fourier integrals 

J] 5 (El{z)A»(z)) = I Vb s exp \i f dz b s (z)E s J(z)A»(z)\ (s = 1,2), 
xes J ^ JzeS J 

(20) 

where b s (z) is a field 'living' on the surface S. It corresponds to the Lagrange 
multiplier in the canonical quantization approach. According to ( |i4|) we turn to 
the calculation of the generating functional of the free Greens functions 

Z^(J,Tj,ij) = C JdAV^V^ Vb exp {i[S^\A) 

+£ (0) (^,V>) + I dz b s (z)E°\z)A»(z) 
JzeS 

+ Jdx (A^x)J^(x) +^(x)ri(x) + fj{x)ip{x)))} (21) 

with SW(A) = lfdx A„{x)K^A v {x), S^{$^) = Jdxtf(x) (id-m)ip(x) and 
K^ v = g^ v d 2 — (1 — l/a)d^d u (hereafter, the summation over double latin indices 
from 1 to 2 is understood). Since the boundary conditions are gauge-invariant, 
they do not interfere with the integration over the gauge group. Therefore, the 
quantization of the gauge field in covariant gauge with gauge parameter a pro- 
ceeds in a standard manner, e.g., with the help of the Faddeev-Popov procedure. 
In order to diagonalize the quadratic form in the exponential of ([H]) we substitute 

A^x) - A^x) - [dy D^(x - y)J v {y) - f dz D/(x - z)E s u (z)b s (z) 
J JzeS 

where 

D» v (x -y) = (g, v - (1 - a)d Xfi d x Jdl) D(x - y) , (22) 

d 4 £; Q ik{x-y) 

W) 

are the photon propagator in covariant gauge and the causal propagator, respec- 
tively. The result is 



r d 4 k e iHx-y) 

= (E>0) (23) 



S {0) {A) + [ dz b s {z)E s ^{z)A^{z) + fdx A^(x)J^(x) 
JzeS J 

l -Jdx A^x)K^A v {x) - l -fdxfdy J»(x)D^(x - y)r{y) 

- ( dx ( dz J fX (x)D u u (x - z)E s u (z)b s (z) 
J JzeS 



If dz f dz' b s (z)K st (z,z')b t (z') 
l JzeS Jz'eS 
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with 



K st (z, z>) = E s Uz)D^(z - z')Ei(z'), z,z'eS. 



- - , -~ ■ (24) 

The newly defined object K st (z,z') is in fact the projection of the propagator 
D fiu (z — z') on the surface S and, with respect to the Lorentz indices, into the 
tangential subspace spanned by the polarization vectors E^(x), s = 1, 2. Further 

~~ 1 st 

we need to define the inversion K (z, z') of this operation 



dz" K- lSt (z, z")K v \z", z') = 6 s (z - z')5. 



JI\ T^t'tl Jl J\ 



St ) 



(25) 



where 5s(z — z') is the delta function with respect to the integration over the 
surface S, 5 st is the usual Kronecker symbol. The quadratic form in fl21~|) is 
finally diagonalized by the substitution 

b s (z) -> b s {z) - J s dz'J dx K- 1 (z,z')El\z')D^(z'-x)r(x) 

and we obtain 

S {0 \A) + J^dz b s (z)E;(z)A fl (z) + J dx A„(x)J"(x) 



dx A^x)K^ A v {x) 



2Js 



dz / dz' b s (z)K s \z,z')b t (z') 



-\jdxj dy J»{x) s D flu (x,y)r(y) 

with the new photon propagator 

s D tlv (x,y) = D^x -y) - D^ u (x,y) 

= D^{x-y) (26) 

-j s dz j s dz> D/(x-z)E;,( Z )K- 1 (z,z')Et,\z')D^ v (z'-y). 
Now, the integration over all fields can be carried out and we arrive at 
Z {0) (J,7J,r]) = C(detK)'^ (dettf)"* 



x exp 



dx dy J»{x) s D, u (x,y)r(y) 



+ - dx dy ip(x)S(x - y)^{y) 



(27) 



where S(x — y) is the usual electron propagator. 

The generating functional ( p7|) together with relation (14) defines the pertur- 
bation theory of QED in the conventional way. Apart from the factor (det K^j 2 
due to the integration over the field b s (z) (z G S), the boundary conditions 
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manifest themselves in the modified photon propagator (p6[). Our calculation of 
the radiative corrections relies on this representation of the boundary dependent 
propagator. As a consequence of the usual language of Feynmann graphs 
can be adopted. It may be observed with the help of fll5f) that the boundary 
dependent part of the propagator D^ u (x,y) and its part depending on the gauge 
parameter a are orthogonal. 

Equation (p7|) was first derived in ]|] in a more complicated way for the case 
of two parallel planes. Here it is derived for a general surface S. However, the 
main problem in a practical calulation consists in finding the inverse of K st (z, z'). 

The representation (|26|) allows for a simple interpretation of the boundary 
dependent part D^ u (x, y) of the photon propagator. It can be viewed as describing 
the free propagation of a photon from the point x to a point z G S on the surface 
S, a subsequent propagation on S from z to z' mediated by K~ x {z, z') which is 
non-local in z and z', and another free propagation from z' to y. In this way, the 
interaction of a free photon with the boundary S is accounted for, as symbolically 
indicated in fig. 1. 

■n 




Figure 1: Illustration of the boundary dependent part D^i^x.y) of the photon 
propagator. 



3 Radiative Corrections 

Using the above derived representation (|27|) we find the effective action by stan- 
dard methods in the form 

r(A, = -THogif + -Trlogi? + S (0) (A) + S (0) ($,^) -z^{lPI graphs} . 

2 2 (28) 

The Feynman rules are obtained from the usual ones (i.e. without boundary 
conditions) by replacing the photon propagator by its modified analogue fl26|) . 
Since the boundary conditions are static, the effective action is proportional to 
the total time T and the ground state energy is given by 

E = -^T(A^^) l _ . (29) 

If there are translationally invariant directions (e.g., parallel to the plates), the 
relevant physical quantity is the energy density and we have to divide by the 
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corresponding volume too. Graphically, the ground state energy is given by 
1PI graphs without external legs. The only two-loop diagram contributing to 

the effective action is 4 > /■ So up to two loop order we find the following 
expression for the ground state energy 
E = E^+E^ = -^TrlogK-^Trlogi^ (30) 

+ ifl dx l dy \- Du ^ x ~ y) ~ Du ^ y) \ u>lu ( y ~ x) ' 

where 

IP" (a; -y) = -ie 2 TrY S c (x - y)"f v S c {y - x) (31) 

is the known polarization tensor. The term proportional to Tr log K and the con- 
tribution of the free propagator D /JiV (x — y) just constitute the two loop ground 
state energy of QED without boundary conditions. These terms do not depend 
on the geometry. According to the renormalization scheme outlined above, their 
removal can be absorbed in a redefinition of the cosmological constant. Therefore, 
they will be omitted in the following. It may be observed that the remaining con- 
tributions are independent of the gauge parameter a, i.e. the gauge dependence 
of the ground state energy is completely contained in its free part. Let us note 



that the representation fl30|) which was first used in [19] is much simpler than 
the one previously applied in || where the ground state energy was expressed in 
terms of the energy-momentum tensor. 

Next we insert the definiton (p6|) of D Ufl (x,y) into the expression ( |30"D for the 
radiative correction and obtain (still for a general surface S) 

E$ = -^fjdzJ dz' K-' st (z, z>) W°(z', z) , (32) 
s s 

where the abbreviation 
W s (z',z) = JdxJ di/^C^i^C^n^Cy-a;)!?/^)^^),^ (33) 

has been introduced. Using the transversality of H^ v (x), 

IV (x) = (g^dl - d x ^d Xv )Il(x 2 ) , 
this expression simplifies to 

tl ts (z',z) = Jdxjdy E%')D{z'-y)g^dln{{y-xf)D{x-z)Ei{z\^ s . 

Introducing the Fourier transform of II(a;) 

n(x 2 ) = / e**» fl(g 2 ) (34) 
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we rewrite it in the final form 

tt'V, z) = El\z>) J ^e^g^^Et(z\ zfes , (35) 

which is best suited for the calculation of the radiative correction in particular 
geometries. There is yet another important conclusion to be drawn from the 
representation (|3~5). The scalar part II(g 2 ) of the polarization tensor is known 
to possess a logarithmic divergence which is independent of q (for example, in 
Pauli-Villars regularization it is of the form — (2a/37r) log(M/m) where M is 
the regularizing mass). However, if only a g-independent constant is added to 
IT, i.e. n(g 2 ) — > n(g 2 ) + c, the quantity Tl ts changes according to U ts (z > , z) — > 
Tl ts (z',z) + c K ts (z',z). Then it may be verified with the help of (|25 ) that the 
corresponding change in the ground state energy is a simple constant which is 
independent of the geometry. Consequently, the removal of the divergence of the 
polarization tensor can be interpreted as a renormalization of the cosmological 
constant in complete analogy to the free Minkowski space contributions discussed 
above. As a result only the finite, renormalized part of the polarization tensor 
needs to be taken into account when calculating the boundary dependent part of 



the radiative correction E 



(i) 



4 The ground state energy in plane geometry 

In this section we wish to demonstrate the calculational techniques needed in the 
spherical case on a simple example first. For this purpose we rederive the results 
for the ground state energy and the radiative correction E^ in the geometry 
of two parallel conducting planes. In this case, the surface S consists of two pieces. 
A coordinization of the planes is given by z = {x a , x 3 = aj}, where the subscript 
i = 1, 2 distinguishes the two planes and a = 0, 1, 2 labels the directions parallel 
to them (they are taken perpendicular to the X3-axis intersecting them at x% = 
|di — <i2 1 = L is the distance between them). The polarizations E s (s = 1, 2) are 
given by @ and do not depend on x a or i. Therefore they commute with the 
free photon propagator (p2|). Inserting them into ( p4|) yields the operator K st in 
the form 

K st (z,z') = -5 st D(x-x% xleS . (36) 

We proceed with deriving a special representation of the scalar propagator. It is 
obtained by performing the integration over ks in eq. (p3|) 

H ~ "ST ' (37) 

with T = \Jk$ — k\ — k\ + ie. Substituting x^ = di and x' 3 = dj we get 

K s \z, z>) = -5 st J ^ ± h l3 e <M--<0 (38) 
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where the abbreviation 

hij =^\ ai - aj \ (ij = l,2) (39) 

has been introduced. With (|38|) we have achieved a mode decomposition of 
the operator K st on the surface S. As an advantage of this representation the 
inversion of K st , defined by (p5|) , is now reduced to the algebraic problem of 
inverting the (2x2)-matrix h-ij. The inverse of hij is given by 

and we get 



K- 1S \z, A = S« / ^ T ^ e ika(xa - X ' a) ■ (41) 

After inserting this expression into ( p6|) we find the photon propagator for the 
electromagnetic field in covariant gauge with conductor boundary conditions on 
two parallel planes as proposed in p. 

Next we substitute (|38D into (|30|) and obtain for the ground state energy 
density per unit area 

^--^WwW-'-s*')- (42) 



where we have divided by the volume V\\ of the translational invariant directions 
parallel to the planes. We perform the Wick rotation k — > ik (thereby T — > 



27 = iyj Izq + kj + and calculate the 'Tr log' 

Tr log U^-L/^ = logdet ^-8 st ±-hiA = 2 log (l - e~ 2 ^) - 41og(2 7 ) . 

The term —4 log(27) yields a distance independent contribution and will be 
dropped. In this way the known result for the Casimir energy is reproduced 

e 



Now we turn to the radiative correction. The polarization vectors E s (s = 
1, 2) commute with II and we find 

We substitute this expression into the radiative correction (|32| ) and get 

- v / d4fc r /r 1 e - ifc3 ^- a ») (U) 



(2tt)< 
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with T = \Jk a k a + ie and k 2 = k^ . Using (HOT) the last equation takes the form 



E, 



(i) 



2i 



d A k 



-iTL 



cos k-xL 



fl(k 2 

k 2 



(45) 



(2tt) 4 sinTL 

By means of the trivial relation exp(— iTL) = exp(iTL) — 2isinrL we separate 
again a distance independent contribution that will be omitted. Now we perform 
the Wick rotation and obtain the final result for the radiative correction to the 
ground state energy to order a in the geometry of two parallel conducting planes 



E 



(i) 



d 4 k 



7 



- 7 L 



cos k-iL 



fl(k 2 



(27r) 4 sinh 7 L V~ *"J k 2 ' ^ 

It is interesting to study the radiative correction ( [46]) in the limit X c /L << 1. 
For this purpose we transform the integration path of the fc 3 -integration. Due to 
n(&3 + 7 2 ) there is a cut with branch point k 3 = i\JArri 2 + 7 2 in the upper half 
of the complex /^-plane. The discontinuity of the one loop vacuum polarization 
n(A; 2 ) across the cut, discIl(A; 2 ) = Yl{k 2 + ie) — Ii{k 2 — ie), is well known 



discfi(A; 2 ) 



2i 



-a\ 1 



Am 2 



1 + 



2m 2 
fc 2 " 



(47) 



As shown in fig. 2 we move the integration contour towards the imaginary axis 
so that the cut is enclosed. The result can be written in the form 



E 



(i) 



d 3 k 



7 



disc n 



itL 3 J (27r) 3 sinh7L J k 3 M m iL 2 k 2 + 7 2 



_ - x /4m 2 L 2 fe|+7 : 



(48) 

In the limit mL » 1, the second term in the round brackets is exponentially 



o7 



7 



iy/Am 2 + q' 2 



Figure 2: Integration path in the ^-plane 



suppressed and can be neglected. The desired series in inverse powers of mL is 
now simply achieved by expanding the square root in the denominator 



E, 



(i) 



d 3 k 



7 e" 



- 7 L 



27rmL 4 7 (2tt) 3 sinli7L 



dk : 



1 k 2 



1 + 



\mLk 2 



(49) 



15 



After some elementary integrations and with the help of 

diSC n (S) 3Z7T 



1 



we find the leading order contribution to the radiative correction 

n 2 a / 1 ' 

2560mL 4 + WL 5 , 



& = *£n + o(-L) (51) 



in agreement with || and |19| 



5 The radiative correction in spherical geome- 
try 

This section is devoted to the calculation of the radiative correction to the ground 
state energy for a perfectly conducting spherical shell of radius R. Appropriate 
coordinates z G S are given by z = {zo,r z = R,9,4>} where r z ,9,<f> are the 
usual spherical coordinates. In order to calculate the radiative correction ac- 
cording to the representation fl32|) we need to know K st (z,z') and Tl st (z, z'), i.e. 
the projections of the scalar propagator D(x — y) and its product with the vac- 
uum polarization, (DH)(x — y) = / dx' D(x — x')H(x' — y), with respect to the 
polarization vectors E s of (|T^). 

Due to the rotational symmetry it proves convenient to express these quanti- 
ties in terms of spherical waves. Using the expansion of a plane wave in spherical 
waves 

= ^^i l j l (gr z )Y^(9 q ,<P q )Y lm (e z ,<p z ) , 
i m 

the scalar propagator (|23|) can be written as 

D(x-y)=J2 J ^e^ -^^^,^)^^;^,^)^^,^) (52) 

lm _oo 

with 

di{uj\r x ,r y ) =i|w|j I (|u;|r<)/i{ 1) (|w|r > ) , (53) 

where ji and are the spherical Bessel functions and Yi m (9,(f)) are the sur- 
face harmonics. r<(r>) denotes the smaller (larger) one of the radii r x and r y . 
The corresponding spherical expansion of {DU) (x — y) is obtained from (|52| ) by 
replacing the quantity di(u;r x ,r y ) by Ui(u;r x ,r y ) whose explicit form is 

2 r°° , 2 -f \ N n(^ 2 -g 2 ) 



z r 00 2 

Ui{u]r x ,r y ) = / dqq ji(qr x )ji(qr y ) 

Ti Jo uo z — q z 

1 El^w^^)^^ 1 (54) 



y/r x Y y ' Jim? 2ixq 2 
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with p = \J q 2 — uj 2 . K l+ i are the modified Bessel functions. 

The polarization vector is proportional to the angular momentum operator 

L (it corresponds to the TM-mode). It therefore commutes with D(z — z') and 
(DU)(z — z'). Consequently, K st (z, z') and Yl st (z, z') are diagonal in s, t and their 
11-components are given by 

K u (z,z') = - [ 7T ^- z ' o) yir n {O z A z )di{u-R,R)Y{ in {d zl , <f> z ,) (55) 

and 

Ii l \z, z') = - J2 J e^*°-^Y lm (9 z , 0,)n,(u;; R, R)Y* m (6 z ,, 2 ,) , (56) 



l>l,m_ 



respectively. Due to the projection with the angular momentum operator L the 
contribution with / = is absent. 

The second polarization vector E 2 contains the operator nx L (it is analogous 
but not identical to the TE-mode). This operator does not commute with either 
D(z — z') or (DH)(z — z'). Substituting the polarization vector E 2 of (|T9|) into 
(II) we find 

K 22 (z,z>) = E%)g^D(z-z')E^ 2 (z') 

J do; e Mzo-z> ) 
~ 1(1 + 1)^-1(1 + l)] dl ^ R > R) 

x [l 2 (l + l) 2 Y lm (0 z , <f>z)Y^, 0,0 (57) 
-u 2 R 2 (n z x L z ) ] Y lm (6 z , faWJOj, $ z <){n z , x L z ,) . 

The method of reduced matrix elements which is known from elementary quan- 
tum mechanics |T7| provides a plain way to simplify the last expression. After a 
straightforward computation we arrive at 

K 22 (z,z')= £ J ^e^^YU0 z A z )fi^,R,R)Y;j9 z/ ,^ (58) 

l>l,m_ 00 

with 

jg + g g) - ^ [g + ih^ + 

* ^ = ^ - 1(1 + 1) • (59) 

Again, the corresponding result for Tl 22 (z,z') is obtained from K 22 (z, z') by re- 
placing di(u; R, R) — > Hi(u] R, R). Since K st (z, z') is diagonal in s, t as well as in 



17 



the quantum numbers {u,l,m} its inverse is readily found. We get 

K- lU (z,z') = y - [ —^ e Mzo~z>) Y ( e ^ -Y*(6 z ,,(p z/ ), 

K-^ 2 (z,z') = £ / e^-^YUOzM f , , o jx KnV*,M -(60) 

Inserting these expressions into (|26|) we find the photon propagator in covari- 
ant gauge with conductor boundary conditions on a sphere. We note that it 
has a somewhat unusual form because it seems not to contain radial derivatives 
(corresponding to Neumann boundary conditions). However, such derivatives 
are implicitly present in the terms proportional to di±\ in fi (|59|) . Utilizing the 
recursion relations for the spherical Bessel functions, fi can be cast into the form 

f t {u,;R,R) = ^k ^ _ ^ - i} , (* = M). 

In fact, substituting K st (z, z') with the last expression into fl3~0|) and performing 
the Wick rotation uj — > we obtain 



^^|(2 i + l)/|.o g (l-A ?( 



x)) (61) 



with x = \y\ and A/ = (sjej)' = (^z+i (x)K l+ i (xfj . This mode sum representa- 
tion was first derived in [TO] with the help of the dyadic formalism. We therefore 
reproduce the known result for the Casimir energy of a conducting sphere. 

Now we return to the calculation of the radiative correction. Substituting 
K~ ist (z, z') and T[ st (z,z') into E^ (|32|) and making use of the orthogonality of 
the surface harmonics one finds 



2,^ J 2tt 



1\ 2 



2 



(62) 



. iMo^R) /(/ + i)n - ^ [m I+ i + (/ + i)n,-i] 



dtico; R, R) 1(1 + 1)^ - fgi + (Z + l)d,_J 

where we have introduced a regularizing factor which is removed by taking the 
limit e — > in the end. Next we perform the Wick rotation uj — > ik. Thereby we 
have from (|53|) 

di(ik;R, R) = ^I l+ i_(\k\R)K l+1 (\k\R) (63) 



and from (54) 

Ui(ik; R, R) = % - J ^ disc% 2 ) I l+k (pR)K l+l {pR) (64) 



4m 2 
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with p = \Jq 2 + k 2 . Finally, we substitute di ( |63| ) and 11; 
correction and obtain 



into the radiative 



E, 



(i) 



2ixR 



E 

1=1 



71 



(k 2 + V 2 ) 



dg 2 



discll 



4m 2 /? 2 



R 2 



where the abbreviations v = I + | 



2 ' 



#1 



h{p)K v {p) 



and 



5? 



hi 

2v 



Uk)K v {k) 



(65) 



(66) 



(67) 



have been introduced. 

Evidently, the expression (|65|) for the radiative correction to the ground state 
energy of a conducting sphere contains ultraviolet divergencies. Hence, it needs 
to be renormalized. In order to isolate the ultraviolet divergent terms we employ 
the uniform asymptotic expansion of the Bessel functions for v — > oo and k — > oo 
with fixed | 

k 2 + q 2 V 2 



1 + 



t 



' k 2S 

l + V 2 . 



(68) 



The corresponding asymptotic expansion of the quantities B\ and B 2 reads 



Bi, 



E E x 



(1,2) * 



2i 



t 2 ' 



+ i-^i i-^i ^ m 3 „2n ' 

1 n=0 i.i>0 47 

n<.i,j <3n 



(69) 



where only even powers contribute. The coefficients X^'P are listed in Appendix 



A. 



Now we insert the expansions 



into 



and perform the integration over 



dk ^-l+2i s l+2< 



^r(^ + ^ + j)r(|; 
r(| + i + i) 



r, / 1 • e ~ 1 ■ 6 • ■ 

x 2 f 1 \- + i, ~y~ + * + j, 2 + 1 + 



In the next step we utilize a Mellin-Barnes representation of the hypergeometric 
function and rewrite the last expression as 



dk t e - 1+2 *s 1+2i 



+ ^ 



(70) 



x 



da ra + i + a)T(^ + i+j + a)T(-a) 



2vri 



m + i+j + a) 
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where the integration contour c goes from —ioo to ioo parallel to the imaginary 
axis with -§ < Sic- < 0. The sum over the orbital momentum I can now be 
carried out. The result is a Hurwitz zeta function: 



oo o 

E v Mn - 2 ° = ( H (2n + 2a + e - 2; |) . (71) 
1=1 1 

Here and in the previous steps we have assumed 3?e > 3 so that 9?e+29for+2n— 2 > 
1 is satisfied. Finally, we put (f47|) into the radiative correction fl65|) where we are 
able to perform the g-integration as well 




3y^F (1 -q)r(-q) 
4 r(f-a) 



(72) 



(we have 3Rcx < along the integration contour c). As a result of all these 
manipulations we obtain the representation 



" Svri? n ^r(z + |) 

da T(l + i + a)T( e -^ + i+j + a)T 2 (-a)(l -a) 



2vri r(| + z + i + (7)r(| -a) 

x (4m 2 i? 2 ) a f2n + 2a + e - 2; . (73) 

With this representation of at our disposal, we are in a position to separate 
the ultraviolet divergent contributions. They are identified with those poles of 
the integrand of ([73]) that cross the integration contour when removing the reg- 
ularization, i.e. in the limit e —>■ 0. All such poles are contained in the addends 
n — 0, 1. Namely, the zeta function possesses single poles at a = (3 — e)/2 (n — 0) 
and a — (1 — e)/2 (n = 1), respectively while the second Gamma function in the 
numerator of ( |73|) contributes a simple pole at c = (1 — e)/2 when n = 0. The 
addends with n > 2 only contain poles that remain in the half-plane 3?cr < 
when taking the limit e — ► 0. The residues of the poles passing the integration 
contour for e — > give extra contributions to Eq . These additional contributions 



are 



r-idiv,l 16 3 2 /- I, 

E = -—am it (74) 
9n 



from o" = (3 — e)/2 and 



from a = (1— e)/2. As we have anticipated, the ultraviolet divergent contributions 
are proportional to R 2 and R°. Hence their removal can be interpreted as a 
redefinition of the parameters a and k in the energy of the classical system (0). 
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The fact that the ultraviolet divergent contributions only appear as finite residues 
is a peculiarity of the zeta-functional regularization. 

The renormalized radiative correction to the ground state energy is now given 

by 



E 



(i) 1 



0)-E, 



div,l 



E 



div,2 



(76) 



a 



8nR 



x 



12 12 [ X nij + X ni. 
n=0 ij ■*■ V " 1 2 < 

da_ £(| + i + a)T(-\ + % + j + a)Y 2 { 
2vri V{i + j + a)V{\-a) 



-(2) 



1 



-a)(l 



x 



(Am 2 R 2 )° ' (h (2n + 2a-2; 



(77) 



At the first glance Eq (|77|) looks just like E'q (|73|) with e = 0. The difference 
is that the poles at o = (3 — e)/2 and <r = (1 — e)/2 are now located on the 
right of the integration contour. The structure of the asymptotic series (|77|) is 
well known. Shifting the integration contour to the left yields an expansion of 
Eq in inverse powers of mR (likewise, an expansion for mR << 1 is obtained 
by moving it to the right). 

Having separated the ultraviolet divergencies with the help of representation 
(|73f ) we return to eq. fl65l). When inserting (|65D into (f76|) the limit e — > 
cannot directly be taken. Therefore we subtract the first two terms (n = 0, 1) of 
the asymptotic expansion fl69|). 



-Bas — , E/ E/ (^nij X nij) 



,2/ 







n=0 i, J 

n<i,j <3n 



2n 



from the integrand and add them again. In this way we split Eq 

,(i) rcn 



(i) r 



Ek 



E f + E a 



(78) 

according to 
(79) 



into 



E f 



and 



E n 



2txR 



tf A i 



V 



2ttR 



E 

i=i 



2tt (k 2 + v 2 ) 



dk v 
2tt (k 2 + v 2 ) 



^discflf-^) (B 1 + B 2 -B as ) (80) 



4m 2 R? 



dg 2 



4m 2 R? 



discn N=- s, 



# 2 



div,2 




51) 
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In the first contribution, Ef, the limit e — > can be carried out under the signs 
of the integral and the sum because they converge now. This form is suited for a 
numerical evaluation, if necessary. In the second contribution, E as we pass again 
to the Mellin-Barnes representation and observe that E as is identical with the 
first two addends (n = 0, 1) of the asymptotic series in (|77|) . 

Representation (|79|) is the final result for the renormalized radiative correction 
of order a to the ground state energy of a conducting sphere. 

Next we turn to the interesting physical situation where the radius R of the 
sphere is large compared to the Compton wavelength of the electron, i.e. mR » 
1. In other words, we are interested in the leading order term when expanding 
Eq in inverse powers of mR. The representation ([T7j) serves us as a starting 
point. As already mentioned above, the expansion in inverse powers of mR is 
achieved by shifting the integration contour to the left. The leading term of the 
desired expansion is given by the residue at a = —1/2. A table of all addends 
contributing a pole at a = —1/2 is displayed in Appendix B. We observe that the 
addends with n = 0, 1, 2 contain also double poles which lead to the appearance 
of logarithmic contributions whereas the terms with n > 2 show only simple 
poles. It is therefore technically convenient to calculate the contributions with 
n < 2 and n > 2 separately. Accordingly, we split the leading order term into 
two pieces, Eq = E n < 2 + E n>2 + 0(m~ 2 R~ 3 ). The residues contributing to E n < 2 
are calculated by standard methods and we find 

a \og(mR) a f 1447573 31 

n - 2 _ 420vr mR 2 87rmR 2 1 16934400 + 5376 
393 . _ 3 u# , 3, 3 u , , 3 



a 



mR 2 



{7.5788 • 10" 4 \og(mR) + 6.4735 • W 3 }. (82) 



The contribution E n>2 can only be calculated numerically. For this purpose, we 
cast it into the more convenient form 



where we have introduced the notation 



Bl = I, T nL (84) 



2v I v {k)K v {k) 

and 



5° = ^ ^1 . ( 8 5) 



t- 2 l 






1 


2uI u {k)K u {k) + k 2 | 




r 4+i(A;)^ + i(A;) + - 

2 


^I v - X (k)K v - X {k)} 
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Expression ( p3[ ) can also directly be obtained from (p0|). For this purpose, the 
term with n = 2 of the asymptotic expansion (|69|) has to be subtracted from 
the integrand in addition to B as . Then, the leading term of the expansion of the 
integrand in inverse powers of q 2 yields E n>2 . 

The remaining task is the numerical evaluation of E n> 2. The result is 

E n>2 = ™ 9.8230 • 10" 6 , (86) 

mR 2 

where contributions up to Z = 50 have been taken into account (the sum over Z 
converges as J2i u ~ 3 , {y = 1 + 1/2)). 

So we finally obtain the radiative correction in leading order l/(mR) 

E^> = 5_ (7.5788 ■ 10" 4 log mR + 6.4833 ■ lO" 3 ) + O ( —!^) • (87) 

mR 2 v > \m 2 R i J 



6 Summary 

In the present paper we have calculated the radiative correction to the ground 
state energy of a perfectly conducting spherical shell. Being proportional to the 
fine structure constant a, it is small compared to the Casimir energy itself. In a 
realistic physical situation where the radius R of the sphere is much larger than 
the Compton wavelenght A c of the electron, the radiative correction is propor- 
tional to the ratio A c / R. In contrast to the case of two parallel conducting planes, 
the radiative correction for a sphere exhibits a logarithmic dependence on X c /R. 

The calculations have been performed in the framework of general covariant 
perturbation theory with the boundary conditions incorporated as constraints. 
A brief review of quantization with boundary conditions in this formalism has 
been supplied. As a valuable advantage of this approach, the ghost degrees of 
freedom are not affected by the boundary conditions so that they need not be 
taken into account. The use of the formalism has been demonstrated by re- 
calculating the Casimir energy including radiative correction for two perfectly 
conducting planes. The Casimir energy of a conducting sphere could be repro- 
duced as well. As expected from the known heat kernel expansion the removal 
of the geometry-dependent ultraviolet divergencies can be absorbed in a redefi- 
nition of the parameters (|7|) of the external system represented by the surface S. 
In the calculation of the vacuum graph to order a the electron mass m and the 
charge e need not be renormalized. 

A thorough investigation of the renormalization of gauge theories in the pres- 
ence of boundary conditions beyond the vacuum structure as well as the anal- 
ogous calculations within the MIT bag model seems to be an interesting future 
perspective. 

The authors thank K. Kirsten and D. Vassilevich for helpful discussions. 
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A The coefficients X 



(1,2) 
nij 



The coefficients X^f in formula (|69|) are: 



(1,2) 
000 



X 



(1) 

1*J 



/ 

i 



1 3 
8 4 



-I o 



8 * 




o o y 



(2) 
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3 
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, 5 
\ 8 











/ 



and 



X. 



(i) 

2ij 



^(2) 
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32 
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32 
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32 
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32 
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21 
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57 


21 
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32 


4 


32 


16 


125 


545 


15 


35 











32 


64 


32 


61 
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16 

















1155 
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B The pole structure of E { 



(i) 
o 



ren 



in ([77]) 



Here we list all addends of ( |77|) which contribute a pole at o = — | when shifting 
the integration path to the left. They result from the zeta function and the first 
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two Gamma functions in the numerator of ([771): 

r(i) r(i + j-l) CH(2n-3;3/2) 



71 


z 




J 




1 [i 















X 


1 







1 




X 


1 







2 




X 


1 







3 




X 


1 


1 











2 





2 




.6 


X 


2 


2 







.4 




2 


4 







.2 




2 


6 











3 . . . oo 





n . 




3n 


X 



x (double pole) 

x (double pole) 



x 



x (double poles) 

x 
x 

X 



The residues can be calculated by standard methods. The double poles yield the 
logarithm log mR. 
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